Landau damping and the echo effect in a confined Bose-Einstein condensate 



00 
On 



A. B. Kuklov 

Department of Applied Sciences, The College of Staten Island, CUNY, Staten Island, NY 10314 

Low energy collective mode of a confined Bose-Einstein condensate should demonstrate the echo 
effect in the regime of Landau damping. This echo is a signature of reversible nature of Landau 
damping. General expression for the echo profile is derived in the limit of small amplitudes of the 
external pulses. Several universal features of the echo are found. The existence of echo in other cases 
of reversible damping - Fano effect and Caldeira-Leggett model - is emphasized. It is suggested to 
test reversible nature of the damping in the atomic traps by conducting the echo experiment. 
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Recent achivements in trapping and cooling alkaline 
gases have raised strong interest to various aspects of 
many body collective behavior. In many respects, prop- 
erties of a confined atomic cloud forming Bose-Einstein 
condensate are unique. For example, low energy col- 
lective modes of Bose-Einstein condensate confined at 
nearly absolute zero temperature demonstrate damping 
whose rate is several times larger than that in the thermal 
cloud §. Now it is commonly accepted that the damp- 
ing in a confined condensate should be Landau damping 
(LD) typical when collisions between quasiparticles are 
rare events It has been especially emphasized in 

Ref. H that the LD occurs in anisotropic traps where the 
spectrum exhibits random character. 

Landau damping (LD) occurs due to a resonance inter- 
action between a collective mode and quasiparticle exci- 
tations. This damping does not lead to a thermalization 
of the quasiparticle distribution. Therefore, it is essen- 
tially a reversible phenomenon of dephasing of the collec- 
tive mode due to mixing with the quasiparticle contin- 
uum in a sense of the Fano effect and Caldeira-Leggett 
model §]. 

A main test of the reversibility of the damping is the 
echo effect. In the conventional phonon echo model (see 
in Ref. [Qj), a phonon mode is characterized by an inho- 
mogeneous broadening. As a result, the amplitude of the 
oscillations initiated at the time moment t = decays in 
time. Second short pulse imposed on this mode at t = r 
creates a partial time reversal of the dynamics causing 
the echo at the time moment t = It. This type of echo 
is called ei The phonon echo is a single mode effect 
whose essential ingridient is a non-linearity of the mode 

The case under consideration is the multimode effect. 
The essence of the echo effect considered here can be 
briefly outlined as follows: externally excited collective 
mode interacts resonantly with other degrees of freedom 
- quasiparticles constituting a thermal component. As a 
result, the amplitude of the mode decays in time. A sec- 
ond external pulse partly reverses in time the evolution 
of the whole system. Hence, a part of the energy returns 
back to the collective mode at the time moment which 
approximately equals to twice the time interval between 
the two external pulses. In accordance with the classifi- 



cation ||, this type of echo can be called e-2- As it will be 
shown below, a crucial element for this echo is a paramet- 
rical excitation of the thermal component. In its nature 
the echo mechanism described above is close to the echo 



in classical uniform collisionless plasma |10|, where it oc- 
curs at the wave vector equal to the difference between 
the wave vectors of the imposed pulses (h]] . Accordingly, 
the delay between the echo and the second pulse can be 
controlled externally. The amplitude of the plasma echo 
p0{ turns out to be of the first order with respect to the 
amplitudes of the both pulses. These properties are due 
to two circumstances: 1) the uniformity of the plasma 
giving rise to the momentum conservation, and 2) the 
quadratic non-linearity of the kinetic equation Jio| . In a 
non-uniform confined condensate, the condition 1) does 
not exist. Therefore, the plasma type echo |l0| should 
not be a universal feature of such a condensate. 

In Ref. jllj the e2-echo effect has been predicted for 
isotropic condensate, where no LD is expected to occur 
|| and where a major mechanism of damping should be 
thermal dephasing pl] |. In this paper the echo in the 
presence of the LD in a confined condensate will be con- 
sidered. It will be shown that the e2-echo occurs in this 
case. Its amplitude turns out to be linear in the initial 
amplitude of the collective mode and is quadratic in the 
amplitude of the second pulse. The time profile of the 
echo response demonstrates a specific double peak struc- 
ture. 

In order to describe the LD and the echo effect in the 
case of a many boson system demonstrating the phe- 
nomenon of Bose-Eisntein condensation in a confined ge- 
ometry, the standard form 



H = J dr&[Hi + 



H 



2m 



Ve, 



(1) 
(2) 



of the Hamiltonian is employed. Here the Bose operators 
^, v&t obey the usual Bose commutation rule; n is the 
chemical potential; V ex denotes the trapping potential; 
g stands for the interaction constant; and m is atomic 
mass. In what follows the atomic units (J% = 1) will be 
employed. 
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Below we will closely follow the approach employed 
in Ref. ||. Accordingly, two types of averaging (...) and 
(■■■)eq are introduced. The first one is performed over the 
initial state, while the second is the equilibrium thermal 
averaging. Thus the Bose field 'J, the excited condensate 
wave function <E> and the equilibrium condensate wave 
function <£>o can be represented as 



(3) 



Here "J stands for the non-condensate part. In order to 
simplify the following consideration, the Hartree approx- 
imation for the thermal component described by \& will 
be employed. As it is well known, this approximation 
is valid at high temperatures [Q. An extension to gen- 
eral case will be considered elsewhere. Consequently, we 
ignore the anomalous mean (vP^/) and retain only the 
quantity n = ($^(x)\I'(x)) in the equation for the con- 
densate wave function <f>. 

Equations for ^ have be obtained in Ref. |§. Here 
a Hartree limit of these will be utilized. For this purpose 
the Heisenberg equation ih4> = [\&, H] is averaged over 
(...) and the three operator mean (^T^r^) i s factorized 
after the substitution of the representation (|J) has been 
made [§. Then, the factorization *t\jr\|r _> 2(* t *)^ + 
(\]/Yl/)»]>t H i s made in the equation for ^. Finally, in the 
Hartree limit these equations are 



ih<& = [Hi + g\$\ 2 + 2gh}§, 
iW = [Hi + 2. 9 (|$| 2 + n)]#. 



(4) 
(5) 

Eqs.(||), (||) together with the definition of n deter- 
mine completely the mean field dynamics of the con- 
densate and the thermal component in the Hartree limit 
(T >> n) in the collisionlcss regime. In this regime the 
only mechanism of damping is the LD which is rather a 
reversible dephasing of collective oscillations than their 
irreversible thermalization. 

It is enough to consider linearized dynamics of the 
quantities <£' = $ — <I>o, n' = h — n°, describing small de- 
viations from the equilibrium, where the notation h = 
{^^) eq is employed || . In the Bogolubov representation 
for the condensate part and the Hartree approximation 
for the thermal component 



(u n a n + w*a*), 



(6) 
(7) 



a n , a* are classical numbers and the operators b n and 
b\ destroys and creates, respectively, thermal particle on 
the level with the energy e n . The values u n ,v n form the 
eigenvector of the Bogolubov system 



H[u n + g$lv 
H[v n -I 
H[ = H x + 2« ? (|$ c | 2 + n u ) 



e n v n = H[v n + g®Q 2 u n , 

2 i ~0\ 



(8) 
(9) 
(10) 



and obey the orthogonality condition 



dx[u^(x)u„(x) - w* l (x)w„(x)] = Sr, 



(11) 



For the thermal part in the Hartree approximation, one 
should set the "hole" part ~ v n to zero. Then employing 
Eqs.(f|), (11), one obtains 



+ 2g A rnk l f, 

where the notations 



(12) 



/./ 



A m ki= / dx$ «„ + v* m )u* k u u (13) 



fki = (t>lbi}-fki, C = (blb k ) eq = f^S, 



kl- 



(14) 



are introduced, with — l/(exp(£fe) — 1) being the 
equilibrium pop ulation at the /c-th level. We employ 
(14|) and obtain the linearized equation 
p ) as 



uibrium pop i 
Eqs.(|,©, @, 
for f kl (see Ref. [ 



ifkl = unfki + 2g(/fc 0) - // 0) ) ^2(A* mkl a m + A m i k a* m ), 

m 

(15) 

where aj kt = e\ - e k . 

Let us assume that an external resonant drive, imposed 
on the system at times —00 < t < 0, has prepared a state 
with some ai(0) ^ and the rest a n (0) = 0, n ^ 1. We 
also assume that the thermal component is not affected 
by this resonant drive, so that /w(0) = 0. The evolution 
between the time moments t = and some t = r is 
characterized by Landau damping due to the terms ~ f k i 
in Eq.(12), so that the amplitude of the oscillations set by 
the drive at t — will decay exponentially at later times. 
The evolution of the excited mode cii (t) can be obtained 
from Eqs.([l2|),(|l5|). In the lowest order one should ignore 
all the low energy modes but ai(t). Then the rate of the 
exponential damping follows as M jl — 7l(£i) where 



(0) 



/, (0) )«(« 



£1 +e k )\A r< 



kl 



kl 



(16) 



It is worth noting that Eqs.p2), (fLq) can be inter- 
preted in a sense of the Fano effect Jjfl . Indeed, the state 
£1 can be considered as a discrete state mixed with a 
quasi-continuum of the pair excitations characterized the 
spectrum u k i. As a result, the discrete state autoionizes 
with the mean life time @- 

Now let us consider the implications of imposing a 
short external pulse at the time moment t — r since 
the state with ai(0) ^ was created. Let us assume 
that r > I/7L so that the solution ai(t) has apparently 
dissipated at the moment of the pulse. This pulse can 
be introduced by means of varying the external potential 
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V ex , so that V ex is replaced by V ex + 8V ex (t) in H x (g) 
with 

SV ex = V 2 (r)5(t- T ), (17) 

where V2{r) denotes some time independent coordinate 
function describing strength of the external pulse. 

Below it will be shown that this pulse will cause the 
echo - a stimulated revival of the mode a\(t) in the 
form of a double peak centered close to the time moment 
t = 2r. Note that, in contrast with a standard plasma 
echo pfj| which essentially requires a non-linearity, the 
echo in the considered case is expected to exist in the 
linear approximation given by Eqs.()l2]),([i"5"|), if the para- 
metrical excitation of the thermal component by the sec- 
ond pulse is taken into account. This excitation produces 
a change of fkiif) while t crosses the moment t = r, so 
that fki (t + e) acquires an admixture of its complex con- 
jugate /fci(r)* = /ifc(r). This fact constitutes a partial 
time inversion of the original evolution of the solution 
ai(t). Considering the solution of Eq.(^), where Hi has 
been modified by (|l7j), in the limit V2(r) — > and ig- 
noring the non-linearity of Eq.(||) during the action of 
the 5-pulse @, one finds #(r + e) = S[V 2 ]V(t). Here 
S[V 2 ] wl - iV 2 stands for the evolution operator trans- 
forming the solution from the time just before the pulse 
to the time just after the pulse. Employing the orthogo- 
nality condition ([ll]) and Eq.(^) in the Hartree approx- 
imation, one can find b m (r + e) = ^ n 5 l mn 6„(r), and 
finally Eq.Q yields 



ks 



^mk^nsfks{ T )j 



S mn = \ d*U* S[V 2 ]u n . 



(18) 



(19) 



The term responsible for the echo corresponds to k = 
ri, s = m in the sum (|l8|). Therefore, in the lowest order 
with respect to V 2 one finds 



fmn(T + e) = \(V 2 ) mn \ 2 f nm (r) + 

(V 2 )mn = / dKU* m V 2 U n . 



(20) 
(21) 



where £' represents terms which do not lead to the echo 
and which will be omitted in the following analysis. 

It is convenient to introduce a quantity A(t) = ax(t) + 
a*(t). Given the symmetry of the original Hamiltonian 
with respect to time inversion and the representation Jl3] ) 
one finds A* mn = A Vmn = Ai nm . Then Eq.(|l2[) yields 
A(t) = —iei(ai(t) — a*(i)). The echo time profile can be 
obtained by finding the Laplace transforms of (|l2|) , Ji~5|) 
on the time interval t > r and expressing them in terms 
of fmn(T+e) andA(r+e), A(r+e). The echo term resides 
in the part containing fmn(i~+e). Thus, the contributions 
~ A(t + e), A{t + e) can be omitted. Then, / mn (r + e) is 
expressed by means of Eq.(^0|). After that the quantity 
fmn{ T ) can be obtained as the inverse transform on the 



time interval < t < oo with f mn (0) = and A(0) ^ 
0, A(0) ^ as though no pulse at t = t is present Q. 
In the solution A(t) only that part A ec (t) which contains 
the e2-echo should be retained. Finally, one finds the 
solution A ec (t) for t > r as 



A ec (t) 



8ef 



de 



(cj -e)(uj' + e) 



1l {lu) 1l {oj')[iA{0)+uu'A{Q)] 



- e\Y + (2e m (u;)) 2 ][(a;' 2 - e 2 ) 2 + (2e llL {uji) 



(22) 



where the notation 



1(e) 



WE 



\(V 2 ) nm \ 2 Af mn (f^ ~ fPMe -s n + e m ) 



(23) 



is employed. In fact, all the information about the pulse 
( p7| ) is contained in the quantity 7 (p3|). 

The double singularity e){u/+e) in Eq.(f22|) gives 

rise to the term — ir 2 5(uj +lj')S(s —u>) plus regular terms, 
in accordance with Ref . jjj . Omitting these regular terms 
and performing trivial integrations over e and a/, one 
arrives at Eq. (|22|) rewritten as 



7T 



due 



-iu(t-2r) 



7H 7L H 2 [ii(0)-q;A( 0)] 

[ (w 2_ £ 2 )2 + (2£i7lM) 



(24) 



212 ' 



This equation is the main result of the present work. 
First, it indicates that the echo is of the e 2 type ||] due 
to the exponent exp(iw(i — 2r)). Second, the integrand 
Eq.(p4|) has complex poles of the second order. This de- 
termines a specific double peak structure of the echo re- 
sponse. Indeed, setting A(0) = 0, Jjiu) = 7 l, 7l(w) = 
7 l(ei) = 7 , one finds the integral (|24|) as 

A ec (t) = 2A ^ (t - 2r)e-^ |4 - 2T| sin( £l (t - 2r)), t > r. 



£1 



(25) 



Note that the echo becomes maximal at \t — 2r| ~ 1/ 7 l 
and reaches the value A max ~ .4(0)7/71, which does not 
depend explicitly on the interaction constant g, as a com- 
parison of Eqs.(p^), (16) indicates. On the other hand, 
A max is proportional to the first power of the initial am- 
plitude of the collective mode A(t) at t = and is of the 
second order with respect to the amplitude of the second 
pulse V 2 — *■ 0. These features are universal and do not 
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depend on details of the matrix elements. The only re- 
quirement is that the sum (|23| ) is well defined at e = £1. 
In anisotropic traps the spectrum of the pair excitations 
can be considered as a continuum p] . Therefore, this sum 
can be replaced by integration which yields a finite value 
of 7. In contrast, a temperature dependence of A max is 
sensitive to details of the matrix elements in Eqs.(p3|), 
©• 

In the presence of irreversible dissipation determined 
by the higher correlators and characterized by some 
rate 7(j rr ) the echo amplitude is to be exponentially 
suppressed. It is straightforward to realize that the 
maximum echo amplitude acquires an extra factor 
exp(— 2r7(i rr )). Therefore, in order to observe a distinct 
echo, the condition I/7L < t < l/~/u rr \ should be ful- 
filled. 

It is worth noting that the structure of the echo re- 
sponse represented by Eqs.(j2^), (|25|) is typical for other 
models of reversible damping as well. Let us discuss this 
for the Caldeira-Leggett model M. In this case some 
oscillator Q interacts with a quasi-continuum of the os- 
cillators Xj. The e2-echo discussed above occurs as long 
as a parametrical excitation of the quasi-continuum is 
imposed. Such an excitation can be introduced by the 
time-dependent energy term H'{t) — <Wf (t)X?/2, 
where 8ujf(t) describes the time dependent part of the 
frequency of the i-th oscillator. If at the moment t — 
some perturbation Q(0) ^ 0, Q(0) ^ was created, it 
will decay due to interaction with the bath oscillators Xi 
||. However, if the parametrical drive of the bath fre- 
quencies is imposed in the form Sujf(t) — kiS(t — r), the 
amplitude Q(t) will demonstrate a stimulated revival - 
the e2-echo. The profile of this echo is exactly given by 
Eqs.(p^), (p5|), where 7 ( |23| ) is replaced by a structure 
which turns out to be of the first order with respect to 
hi. Accordingly, the echo amplitude becomes of the first 
order with respect to the amplitude of the external drive. 
A detailed discussion of the echo effect in the Caldeira- 
Leggett model will be given elsewhere. 

In summary, it is shown that the e2-echo can be ob- 
served in a confined Bose-Einstein condensate in the 
regime of collisionless damping. It has a specific two- 
peak structure. The echo effect described above could be 
observed in the traps employed in Refs. 0- The e2-echo 
is typical for the Caldeira-Leggett model as well. 
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